clear all;close all;clc;
load Real

range = 0.1;
lambda = 0;
tol = 1e-7;
% param = warmStart(L, rank);
param.r = 1;
param.k = 6;
[model, llh, rate, R] = EM_main(nos, range, lambda, tol, param);

est =  model.U*model.V';